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Abstract 

By using the method developed in the paper [G.Pantsulaia, G.Giorgadze, On some 
applications of infinite-dimensional cellular matrices, Georg. Inter. J. Sci. Tech., 

Nova Science Publishers, Volume 3, Issue 1 (2011), 107-129], it is obtained a repre¬ 
sentation in an explicit form of the particular solution of the linear non-homogeneous 
ordinary differential equation of the higher order whose coefficients are real-valued 
simple functions. 
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1. Introduction 

In [4 ] has been obtained a representation in an explicit form of the particular solution of the 
linear non-homogeneous ordinary differential equation of the higher order with real-valued 
coefficients. The aim of the present manuscript is resolve an analogous problem for a linear 
non-homogeneous ordinary differential equation of the higher order when coefficients are 
real-valued simple step functions. 

The paper is organized as follows. 

In Section 2, we consider some auxiliary results obtained in the paper [4]. In Section 
3, it is obtained a representation in an explicit form of the particular solution of the linear 
non-homogeneous ordinary differential equation of the higher order whose coefficients are 
real-valued simple functions. In Section 4 we present mathematical programm in MathLab 
for the graphical solution of the corresponding differential equation. 
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2. Some auxiliary propositions 


For n^N,vje denote by [—Z, Z [ a vector space of all n-times differentiable functions 
on [—Z,Z[ such that a series obtained by Zc-times differentiation term by term of the Fourier 
trigonometric series of 'F pointwise converges to for all G [—Z, Z[ and 0 <k <n. 

Let (A„)o<m< 2 m be a sequence of real numbers, where M is any natural number. For 
each Z: > 1 we put 

<Ji=E(-irA2,(^P', (2.1) 

n=0 ' 

( 2 . 2 ) 

11=0 ' 


Theorem 2.1, fZ@]Z, Theorem 3.1, p.45) For m>l, let us consider an ordinary differential 
equation 


2m jn 

n—0 


dx” 


(2.3) 


where 



(2.4) 


and A„ G Rfor 0 <n< 2m. 

Suppose that Aq 0 and af + (o| / 0/or Z: > 1, where (5k and (Ok are defined by (2.1) 
and (2.2), respectively. 

If (Y,ciffiF 2 ,d 2 , ■ ■■) is such a sequence of real numbers that the series 'Fp, defined 
by 


/t=l 


2A() 


(Ck(5k-dk(Ok\ fkTix\ /Ck(Ok + dkak\ . (k%x\ 

y-iFFF-> “H-) + “"(-)• 


(2.5) 


belongs to the class Z,Z[, then 'Fp is a particular solution of (2.3). 

Theorem 2.2. ( Z@]Z, Theorem 3.2, p.45) For m> I, let us consider an ordinary differential 
equation (2.3), where 

/(x)GC[-Z,Z] (2.6) 

and A„ G Rfor 0 <n< 2ni. 

Suppose that Aq 0 and + (o| 0 for k > \, where (5k and (Ok are defined by 

(2.1) and (2.2), respectively. Let (^,c\,d\,C 2 ,d 2 ,- ■ ■) be Fourier coefficients of f and 
("^, Cl, d \, C2, d2 ; • • •) G .^1. 

Then the series *Fp, defined by 


'i'pix) 


Co yi /Ck(5k dk(Ok\ 



Ck(Ok + dk(5k 
^k + ^k 



(2.7) 


is a particular solution of (2.3). 
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3. A non-homogeneous ordinary differential equation of higher 
order whose coefficients are continuous or real-valued step 
functions 

Let consider a partition of [—/,/[ defined by 

We define a differenfial operator 



for G Notice fhaf L(*L) can be rewriffen as follows 



for G where Ind denotes an indicator function. 

For each 5 G N we define an operator Ls by 

1(2<! -I-1 — S'! d" 

Lsi^) = L Ind^ij2^ i(2s+2-s) j ( £ A„(- - - 

for'F gFD(2'«) [-/,/]. 

Lemma 3.1. For each G have 

L{^>) = YlmLs{^). 

S-^oo 


Theorem 3.2. For m>\, let us consider an ordinary differential equation 



(3.1) 


where 



andAffx) G C[—1,1] for 0 <n< 2m. 

Suppose that Ao(x) = 1 and ct|(x) + (of (.r) / 0/or x G [—l,l[ and k>l, where Ok{x) 
and (Ok{x) are defined by 



(3.3) 
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kll 

i:(-irA2„+i(x)(-)2”+i. 

^—r\ ^ 


Suppose that the following conditions are valid: 

(i) (^ ? ^2 1 ^2 ; • * • ) ^ 

(Hi) There is a constant C > 0 such that 
(>ik{x + h) 

a|(x + /i)+ (of(x + /i) al{x) + oil{x) 

i 

Ok{x + h) _ (5k{x) 

a|(x + /i)+ (o|(x + /i) dl{x) + tiil{x) 

Then the function *Po. defined by 


<C\h\^ 

<C\h\^. 


M, / \ , Y- fCkOk{x)-dk(Ok{x)\ (k%x\ 

'i'oW = ^ + 2^ — 27 ——27^ cos -r + 
2 a^(x) + (o^(x) / \ I J 


Ckt)ik{x)+dk(^k{x)\ . fkttx 
-— sin ^ 


V al{x) + oil{x) J \ I J' 

is a particular solution of (3.1). 

Proof. We put 

^ /n<^r /(2i-5) ;(2.,+2-5) r(x) — + 

i=0 I ^ ’ s I L 2 

^ /CkOk C^^^^^ - dk(i>k C^^^^g /k%x\ 

£ ( (-) + 

I ^2(/(^) + (,2(/(^) J l — JJ ’ 

On the one hand, by using the result of Lemma 3.1 we have 


L(lim '¥s{x)) = lim L('L 5 (y) = lim Lsi'Psix)) = lim f{x) = f{x). 

yoo 0^00 O—>-oo 

On the other hand we have 


I lim -'Lo(x))| = lim |'L 5 (x) -'Lo(x)| = lim ^ Ind.iQs 


s-S) l(2s+2-S) , 


w r?- 





1 CkOk{x)-dk(Ok{x)\ 

fktlx' 

POC 1 _ 


al{x) + Oil{x) ) 

z , 


Cktiik{x) + dkOk{x)\ . 

rk'Kx'y. 


<5l{x) + of,{x) )""" 

w“r )\ 
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S-l 






oo 

sup {£ 

-5) ((2s+2-5')j 

\Ck(^ki ^5 dktHki ^ s OkOkix) dktOkix) 

1 + a- 

■) + dkOk{^ 5 CkOik{x)-\-dkakix)\ 

l{x) + (ol{x) 

} = 

1 ^2(/(2.+ 1-5))^^2(/(2.+ 1-5)) O^x) + trxi{x) 


+ 


S-1 

1^1 


5-s« 


sup 




CTyt( 


/( 2 i+l-S)' 

S 


Okix) 




,2( /(2m-S) ^^^2( /(2.+ l-5) ^ Cr2(;c) + (02(x)' 




CO/t( 


/(2S+1-5) 


)) 


(Sik{x) 


dk{ 


^2( K2m S) )^(o2( /(2^+^l 


2(S2i^)+o,;:(a2£^) of(A-) + t^(A-) 


Ci:( 


(0^( 


S ! ' '^k\ S 
l{2s+l-S) 


'-)) 


(iik{x) 


a2(!e£^) + co2(!P£^) al{x) + oyi{x)' 


+ 


)+ 


}< 


s-i 

lim ^ sup 


4l^C 

{£2(h| + K|)^)< 

k^l ^ 


qi2(~' 

j™-^L(k-t| + |4|) = o. 


k=\ 


□ 


Remark 3.3. Theorem 3.2 is a generalization of Theorem 2.2. Indeed, Theorem 2.2 is a 
simple eonsequenee of Theorem 3.2, when A„(x) = const for 0 < n < 2m, beeause in that 
eases all eonditions of Theorem 3.2 are fulfilled. 

We say that {ak)o<k<s is partition of [—/,/[ if —/ = ao < ai < • • • < a^-i <as = 1. 

We say that a real-valued funetion / on [—/,/[ is simple funetion if there exists a parti¬ 
tion {ak)o<k<s of [—l,l[ and a sequenee of real numbers {Ak)i<k<s sueh that 

fix) = 

k=i 


for X G [—//[. 

We have the following proposition. 

Theorem 3.4. Suppose that (A„(x))o<n<2m is a sequence of real-valued simple step func¬ 
tions on [—/,/[, i.e. for every n {0 <n < 2m) there exists a partition ia^k^)o<k<s„ of [—/,/[ 

in) 

and a sequence of real numbers (A^ )i<k<s„ such that 


^nix) = Y^Af^Ind 




k=\ 
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forx G 

Suppose that Ao(x) does not remain a zero value on [—/,/[ and CT^(^) + (0^(x) 7 ^ 0 for 
X G [—l-,l[ and k>l, where Ok{x) and are defined by (3.3) and (3.4). Suppose also 
that Fourier coefficients of the function f standing in the right side of the equation (3.1) 
satisfy the following condition (‘lf,c\ffi\,C 2 ,d 2 -,---) ^t\- 
Then the function *Po. defined by 




Co 


dktOk(x)'\ / kTlx\ 

+£( )‘°n—J+ 


k=i ' W) + ^ki^) 


2Aq(x) ' V 

( Cktiik{x) + dkOk{x) 


\ I 


sin 




V al{x) + oil{x) J \ I J' 

for X G [—l,l[, satisfies (3.1)-(3.2) at each point of the set 

{-l,l)\ Uo<„< 2 m{al"^, « 2 ”V • • ,«!”-l} 


(3.7) 


Proof. Ifxo G {—fl)\G(\=iJQ<n< 2 m{a^\\a^ 2 ^,--- then by virtue of the openness 

of the G there exists a positive real number r > 0 sueh that (xq — r,xo + r) C G. It is obvious 
that (A„(x) is eonstant on (xq — r,xo + r) for 0 < n < 2m. We set A„ := A„(xo) for 0 < n < 2m. 
For m > 1, let us eonsider an ordinary differential equation 


2m 

14 

n—0 


d" 


(3.8) 


Note that for (3.8) all eonditions of Theorem 2.2 are fulfilled. Henee the series 
defined by 

Co , / CkOk-dkti>k \ tkiix-x 


^p{x) = 


2Ao 


+ I 

k=\ 


V 


ctf + cof 




\ I 


+ 


/ Ck(iik + dk(^k \ . tkllx\ 


V 


Ctf+ cof 


V I J 


(3.9) 


is a partieular solution of (3.8), where <5k and tiik are defined by (2.1) and (2.2), respeefively. 

Nofiee fhaf defined by (3.9) eoineides wifh \|/o defined by (3.7) af all poinf x G 
(xo — r,XQ + r). Similarly, fhe equation (3.8) wifh (3.2) eoineides wifh fhe equation (3.1) 
wifh (3.2) af all poinf x € (xq — r,xo + r). Henee \|/o defined by (3.7) satisfies (3.1)-(3.2) 
af eaeh poinf of fhe sef (xo — r,xo + r), in partieular, af poinf xo. Sinee xo G (—1,1) \ G was 
faken arbifrary, we end fhe proof of Theorem 3.4. 

□ 


4. On a graphical solution of the linear non-homogeneous or¬ 
dinary differential equation of the higher order whose coeffi¬ 
cients are real-valued simple step functions 

Lef eonsider fhe linear non-homogeneous ordinary differential equafion of fhe 22-lh order 

j2 j5 j14 j20 

'T(x) +A 2 (x)^'F(x) +A 5 (x)^'F(x) +Ai 4 (x)^'T(x) +A 2 o(x)^'F(x)+ 
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(4.1) 


where 


Miix) 




*F(x) = 1 +2cos(x), 


A2{x) = —0.001 “0-002 x/?iJ[_;j/2,o[(-^) “0-001 ^^^'^[o,jc/ 2[(-^)“ 0-002 x/?i(i[;j/2,jt[W) 


A5(.r) =0.01 X/?i(i[_;^ _;^/2[(x)-0.01 x/?i(i[_^/2,o[ W +0-002 x/?iJ[o,;;c/2[(x)-0.002 x/?i(i[„/2,;;c[(x), 


Ai 4 (x) = 0.1 X -jt/2[ (.r) - 0.1 X Ind[_T^j 2 fl[ {x) - 0.4 x /?ifi?[o, 7 i/ 2 [ (-^) +0.007 x Ind\j^/2^n[ i^) > 


A 2 o{x) = -0.01 X -jc/2[(-^) +0.01 X /?i(i[_:;t/ 2 ^o[(-^) + 0-002 x Ind[o^n/ 2 [{x) -0.22 x /?i<i[7c/2,7i[(-x), 


A 22 {x) = 0.001 x/?i(i[_;j _;j/ 2 [(-^) “0.001 X /?ir/[_^/ 2 ,o[W + 0-0003 x/?i(ip, 7 i/ 2 [(-^) “0-0003 x /?i(i[^/ 2 , 7 i[(-^)- 
Definition 4.1 We say that g G 7t,7t] \ G) (G := {—7t, —7t/2,0,7t/2,7t}) if 

g{x)=gi{x)y<lnd^_T,_T,l2[{x)+g2{x)xlnd\^_T,l2fi[{x)+g2{x)xlnd[Q^T,l2[{x)+gA{x)y<lnd\j,i2^T^[{x) 

(4.2) 


for some gi,g 2 ,g 3 ,g 4 S FDp2)([-7i,7t[). 

Below we present the program in MathLab which gives the graphical solution of the 
differential equation (4.1) in the class 7t,7t[\G). 


A1 = [0,-0.001,0,0,0.01,0,0,0,0,0,0,0,0,0.1,0,0,0,0,0,-0.01,0,0.001]; 

A2 = [0, -0.002,0,0, -0.01,0,0,0,0,0,0,0,0, -0.1,0,0,0,0,0,0.01,0, -0.001]; 

A3 = [0, -0.001,0,0,0.002,0,0,0,0,0,0,0,0, -0.4,0,0,0,0,0,0.002,0,0.0003]; 

A4 = [0, -0.002,0,0, -0.002,0,0,0,0,0,0,0,0,0.007,0,0,0,0,0, -0.22,0, -0.0003]; 


CO = 2;A10 = 1;A20 = 1;A30 = 1;A40 = 1; 

X = 1 : 20; 

51 = A10; 52 = A20; 53 = A30; 54 = A40; 
for ^ = 1:11 

51 = 51 + (-1) ()t) * A1 (2 * )t) *x. (2 * )t); 

52 = 52 + (-1) (A:) * A2 (2 * A:) * X. (2 * A:); 

53 = 53 + (-1) (A:) * A3 (2 * A:) * X. (2 * A:); 

54 = 54 + (—l)(fc) *A4(2*^) *x.^2*k); 

end 

01 =A1(1);02 = A2(l);03 = A3(l);04 =A4(1); 
for ^ = 1 : 10 
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01 = (91 + (-l))t*Al(2*A:+l)*x.(2*/c+l); 

(92 = (92 + (-l))t*A2(2*A:+l)*x.(2*/c+l); 

03 = (93 +(-l))t*A3(2*A:+1) *x.(2*/c+1); 

04 = 04+(-l)A:*A4(2*A:+l)*x.(2*/c+l); 

end 

xl = {-pi) : {pi/100) : {-pi/2)-, 

yl=C0/{2*Al0)-, 

for « = 1 : 20 

yl = y \ + cos{n*xl). * {C{n) *S\{n) —D{n) * Ol{n))/{Sl{n)2 + 01(?i)2) + 
sin{n*xl). * {C{n) * 01(?i) +D{n) *S\{n))/{S\{n)2 + 01(?i)2); 

end 

x2 = {—pi/2) : {pi/\00) : 0 ; 
y2 = C0/{2*A20y, 
for n = \ :20 

y2 = y2 + cos{n*x2). * {C{n) *S2{n) —D{n) * 02{n))/{S2{n)2 + 02(?i)2) + 
sin{n*x2). * {C{n) * 02{n) +D{n) * S2{n))/{S2{n)2 + 02(?i)2); 

end 

;c3 = 0 : {pi/100) : {pi/2)-, 
y3 = C0/{2*A30); 
for « = 1 : 20 

3^3 = y3 + cos{n*x3). * {C{n) *S3{n) —D{n) * 03{n))/{S3{n)2 + 03(?i)2) + 
sin{n*x3). * {C{n) * 03{n) +D{n) *S3{n))/{S3{n)2 + 03{n)2)-, 

end 

xA = {pi/2) : {pi/100) : pi-, 
y4 = C0/{2*A40)-, 
for « = 1 : 20 

y4 = y4 + cos{n*x4). * {C{n) *S4{n) —D{n) * 04{n))/{S4{n)2 + 04(?i)2) + 
sin{n*x4). * {C{n) * 04(?i)+ D{n) * S4{n))/{S4{n)2 + 04(?i)2); 

end 

for / = 1 : 20 

if01(/)2 + 51(/)2 =0;O2(/)2 + 52(/)2 = 0;O3(/)2 + 53(/)2 = 0;O3(/)2 + 53(/)2 = 


0 ; 


plot(xl , 3 ; 1, x2, 3^2 , x3, j3, x4, 3 ; 4 ) 

else error(’the ordinary differential equation has no solution or has infinitely many 
solutions’ in the class [— 71 , 71 ] \ G)) 

end 
end 

On Figure 1, the graphieal solution of the differential equation (4.1) is presented. 
Remark 4.1 Notiee that for eaeh natural number M > 1, one ean easily modify this 
program in MathLab for obtaining a graphieal solution of the differential equation (3.1)- 
(3.2) in FD^'^^\[—l,l[\G) whose eoeffieients (A„(x))o<n< 2 M are real-valued simple step 
funetions on [—/,/[, / is a trigonometrie polynomial on [—/,/[ and G is the partition of the 
interval [—/,/[ defined by the family (A„(x))o<n< 2 M- 




Figure 1. Graphical solution of the ODE (4.1). 


Remark 4.2 Since each constant c admits the following evident representation 

c = c X Ind[_T^ _T^i 2 \{x) + cx /nd[_^/ 2 , 0 [W + C X + c x (4-3) 

we can use above mentioned program for a solution of the differential equation (2.3)-(2.4) 
with constant coefficients. 

On Figure 2, the graphical solution of the linear non-homogeneous ordinary differential 
equation of the of the second order with real-valued constant coefficients 

'F(;c) - = l/2 + cos(;c), (4.4) 

is presented, which has been obtained by entering in the above mentioned program of the 
following data: 

A1 = [0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

A2 = [0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

A3 = [0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

A4= [0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

C = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

D = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

CO = 1; AlO = 1; A20 = 1; A30 = 1; A40 = 1; 
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Figure 2. Graphical solution of the ODE (4.4). 


Remark 4.1. The approach of Theorem 3.4 used for a solution of (3.1)-(3.2) with real¬ 
valued simple step functions (A„(a:))o<«< 2 m can be used in such a case when the corre¬ 
sponding coefficients are continuous functions on [—/,/[. If we will approximate these 
coefficients by real-valued simple step functions, then it is natural to wait that under some 
’’nice restrictions” on these coefficients the solution obtained by Theorem 3.4, will be a 
’’good approximation” of the corresponding solution. 
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